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We analyze the spatio-temporal patterns of two competing species in 
the presence of two white noise sources: an additive noise acting on the 
interaction parameter and a multiplicative noise which affects directly the 
dynamics of the species densities. We use a coupled map lattice (CML) 
with uniform initial conditions. We find a nonmonotonic behavior both 
of the pattern formation and the density correlation as a function of the 
multiplicative noise intensity. 

PACS numbers: 05.40.-a, 87.23. Cc, 89.75.Kd, 87.23.-n 

1. Introduction 

We present a stochastic model for spatial distribution of two compet- 
ing species. Our theoretical model could be useful to describe biological 
systems, where the presence of fluctuations, such as random variability of 
temperature, can modify strongly the dynamics of an ecosystem We 
focus on the role played by the noise on the transient dynamics of the spatial 
distributions of two competing species belonging to an ecosystem described 
by generalized Lotka-Volterra equations |3j in the presence of multiplicative 
noise. We find nonmonotonic behaviors for the pattern formation and the 
density correlation of the species as a function of the multiplicative noise in- 
tensity. The theoretical results could contribute to select environmental and 
periodical driving forces using the proper space and time scales to develop 
physical models of population dynamics useful to interpret spatial patterns 
in the abundance of the species [U[S]- 

* Presented at the 16 th Marian Smoluchowski Symposium on Statistical Physics: Fun- 
damentals and Applications, Zakopane, Poland, September 6-11, 2003 
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2. The model 

To study the spatial effects due to the presence of noise sources we 
consider a discrete time evolution model, which is the discrete version of 
the Lotka-Volterra equations with diffusive terms, namely a coupled map 
lattice 6 . The time evolution of the spatial distribution for the two species 
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Fig. 1. The bistable potential U(f3) of the interaction parameter [3 (t) . The potential 
U(f3) is centered on (3 = 0.99. The parameters of the potential are h = 6.25 • 10~ 3 , 
7] = 0.05 ; p = -0.01. 



is given by the following equations 
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where xfj and yf j denote respectively the densities of prey 1 and prey 2 in 
the site at the time step n, /i is proportional to the growth rate, D is 
the diffusion constant, indicates the sum over the four nearest neighbors. 
The random terms are modeled by independent Gaussian variables denoted 
by Xfj , Yfj with zero mean and variance unit and (3 n takes into account for 
the interaction between the species. In eqs. Q and ©, a x and a y are the 
intensities of the multiplicative noise which model the interaction between 
the species and the environment. 
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Fig. 2. Time evolution of the interaction parameter (3 for crp — and initial value 
(3(0) = 0.94. The values of the parameters are: 7 = 1.5- 10 _1 , ujq/(2tt) = 10~ 2 . In 
the absence of noise (3(t) oscillates below the critical value (3 C = 1. The transient 
behavior near the initial point (3(0) = 0.94 is due to the choice of cosine function 
in eq. 0) as deterministic periodical driving. 



2.1. Stochastic Resonance 



It is known that for f3 < 1 a coexistence regime takes place, that is both 
species survives, while for (3 > 1 an exclusion regime is established, that 
is one of the two species vanishes after a certain time j7||H]- Coexistence 
and exclusion of one of the two species correspond to stable states of the 
Lotka-Volterra's deterministic model Real ecosystems are immersed in 
a noisy nonstationary environment, so also the interaction parameter is af- 
fected by the noise and some other deterministic periodical driving such as 
the temperature. The change in the competition rate between exclusion 
and coexistence occurs randomly because of the coupling between the lim- 
iting resources and the noisy environment. A random variation of limiting 
resources produces a random competition between the species. The noise 
therefore together with the periodic force determines the crossing from a 
dynamical regime (f3 < 1, coexistence) to the other one (f3 > 1, exclusion). 
To describe this continuous and noisy behaviour of the interaction param- 
eter (3(t) we consider an Ito stochastic differential equation with a bistable 
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potential and a periodical driving force 

d/3(t) dU{[3) . . . (J . 

where U{(5) is the bistable potential (see Fig.0) 

U(0) = h(P - (1 + p)f/^ - 2h(0 - (1 + p)fh\ 



(3) 



(4) 



the periodical driving mimics the climatic temperature oscillations, and 
£p(t) is a white Gaussian noise with < £p(t) >= and < £p(t)£p(t') >= 
apS(t — t'). Since the dynamics of the species strongly depends on the value 
of the interaction parameter, we initially consider the time evolution of (3 
for up = (see Fig.©). We note that in the absence of the additive noise 
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Fig. 3. Time evolution of the interaction parameter f3(t) for op = 2.65 • 10 -3 . The 
values of the parameters and the initial value (3(0) are the same of Fig\^ For this 
level of noise a synchronization appears (Stochastic Resonance) and (3(t) oscillates 
quasi-periodically below and above the critical value f3 c = 1. 



£p(t), (3(t) shows a periodical evolution but its values always remain below 
/3 = 1, i. e. in the coexistence regime. The noise can synchronize with the 
periodical driving force. In this case a Stochastic Resonance (SR) HU] 
effect appears which affects strongly the dynamics of the system |71 ITT1[T^] . 
Therefore we fix the additive noise intensity at the value ap = 2.65 • 10 -3 
corresponding to a competition regime with (3 periodically switching from 
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coexistence to exclusion regions ^2] (see FigEJ). The SR effect in the dy- 
namics of interaction parameter (3 induces SR phenomenon in the dynamics 
of two competing species. This produces noise-induced anticorrelated pe- 
riodic oscillations in the time evolution of the two species densities (see 
ref. [12]. 



2.2. Spatial distributions 

We consider the time evolution of the spatial distribution of the ecosys- 
tem, described by Eqs. (^Q) and (J2J), in the SR dynamical regime. The 
interaction parameter (3 n of eqs.(^Q), © corresponds to the value of contin- 
uous j3(t) of eq.© taken at the step n. So we fix the additive noise at the 
value = 2.65 • 10~ 3 and we vary both the intensities of multiplicative 
noise. In Figs. |1] and [5] we report the spatio-temporal patterns of the two 
species for different values of the multiplicative noise intensity a = a x = o y , 
namely a = 1(T 12 , 10" 8 , 10" 4 , 10" 1 with fj, = 2, D = 0.05, A = 1.5 • 10"*, 
w /(2vr) = 10~ 2 , /?(0) = 0.94 and x° Lj = = 0.5 at all sites (i,j). We 
see that for very low noise intensity (see Figjlji) an average correlation on 
the considered lattice (JV = 100 x 100) between the species is observed. For 
higher noise intensities (see Fig0Ja,[2K) an anticorrelation between the two 
species is observed: the two species tend to occupy different positions. The 
anticorrelation is more evident in FigOJ.. By increasing the multiplicative 
noise the anticorrelation is strongly reduced (see FigJSJa). Further increase 
of the noise (a = 10 +3 ) causes the anticorrelation to disappear and the 
two species densities become uncorrelated. We note also that the average 
size of the patterns increases with the noise intensity non-monotonically: at 
very low noise intensity (a = 10 -14 ) the spatial distribution is almost uni- 
form, by increasing the noise intensity spatial patterns arise (see FigsHJa, 
EK) and a further increase of the noise intensity reduces the average size of 
the patterns (see FigEJa). 



2.3. Spatial correlation 

In order to evaluate the spatial correlation between the two species for 
the noise intensities considered we calculate, at the time step n, the mean 
correlation coefficient < c n > defined on the lattice as 

< c« >= ( 5 ) 

a x°y 



with 
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Fig. 4. Spatial distributions at different times for (a) a = 10 12 and (b) a = 10 . 
The value of the additive noise is fixed at o~p = 2.65 ■ 10~ 3 . The values of the 
parameters are: fi = 2, D = 0.05,. j = 1.5 • 10" 1 , w /(2tt) = 10~ 2 , N = 100 x 100. 
The initial values are j — yfj = 0.5 for all sites and (3(0) = 0.94. 



where x n , s™, y n , s™ are the mean value and the root mean square respec- 
tively of species 1 and species 2, obtained over the whole spatial grid at 
the time step n, cov™ y is the corresponding covariance and N = 100 x 100 
the number of sites which compose the grid. The behaviour of the mean 
correlation coefficient < c n > as a function of the time for different levels 
of the multiplicative noise has been reported in Fig|BJ We observe a non- 
monotonic behaviour of < c n > as a function of the multiplicative noise 
intensity. In fact for low noise intensities a = 10~ 12 , < c n > shows weak 
oscillation around 1, that is strong correlation between the two species. For 
higher levels of the noise a = 10~ 10 , < c 11 > is affected by fluctuations and 
its values vary strongly as a function of the time. A further increase of the 
multiplicative noise, i.e. a = 10 -8 and a = 10~ 4 , determines an oscilla- 
tion of < c n > around a negative value, that is anticorrelation between the 
two species, with the frequency of the periodical forcing. For higher inten- 
sities of the noise a = 10 _1 , the value of the mean correlation coefficient 
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Fig. 5. Spatial distributions at different times for (a) a — 10~ 4 and (b) a = 10 _1 . 
The additive noise intensity, the values of the parameters and the initial conditions 
are the same of Fig^ 



< c n > increases and it vanishes for a = 10 +3 . Finally to show clearly the 
nonmonotonic behaviour of < c n >, we calculate the time average of the 
mean correlation coefficient < c n >t and we report it as a function of the 
multiplicative noise intensity in Fig0 A clear minimum is shown, which 
corresponds to the anticorrelated oscillations shown in the time evolution of 
two competing species in each point of our spatial grid ^21- We note there- 
fore the different role of the two noise sources in the ecosystem dynamics. 
The additive noise determines the conditions of the dynamical regime, the 
multiplicative noise produces a coherent response of the system (see ref.|12p. 
which is responsible for the appearence of anticorrelation behavior in the 
spatial patterns of the species. 



3. Conclusions 

We present a study on the role of the noise in the spatial distributions of 
two interacting species. The main result is that noise plays a very important 
role in population dynamics and cannot be neglected. Noise can have a 
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Fig. 6. Mean correlation coefficient < c n > as a function of the time. For low 
levels of the multiplicative noise (a = 10~ 12 J the species are strongly correlated and 

< c™ > is approximately constant. By increasing the intensity of the multiplicative 
noise (a = lO^ 10 ^ < c™ > shows big fluctuations. A further increase of the noise 
(a = 10~ 8 ; a = 10~ 4 j causes strong anticorrelation between the two species with 

< c" > oscillating at the frequency of the periodical forcing. For very high levels 
of noise, the anticorrelation is reduced (a = lO" 1 ^ and finally it disappears (a = 
10+ 3 ), that is the species are totally uncorrelated. 



constructive role and it is responsible for the enhancement of the response 
of the system to a driving force producing stochastic resonance. By using 
a discrete time evolution model, which is the discrete version of the Lotka- 
Volterra equations with diffusive terms, in the presence of a multiplicative 
noise and with a random interaction parameter, we analyze the temporal 
behaviours of the spatial distributions for an ecosystem consisting of two 
species. The noise induces spatio-temporal behaviors which are absent in 
the deterministic dynamics, i.e. pattern formation with the same periodicity 
of the deterministic force. Moreover appearance of temporal oscillation is 
observed in the correlation coefficient between the two species. We find 
a nonmonotonic behaviour of the time average correlation coefficient as 
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Fig. 7. Time average of the mean correlation coefficient < c n > t as a function of 
the multiplicative noise in semilog scale. 



a function of the multiplicative noise. In fact at low levels of the noise 
intensity the species densities are almost uniform and a strong correlation 
appears in the spatial distributions. By increasing the noise intensity we 
observe pattern formation and anticorrelated behaviour, which exhibit the 
same periodicity of the deterministic driving force and corresponds to the 
minimum of Fig0 For higher values of the noise intensity no patterns or 
correlations appear. Our model could be useful to explain spatio-temporal 
behaviours of populations, whose dynamics is strongly affected by the noise 
and by the environmental physical variables. 
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